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Abstract 

Recent  work  on  the  qualitative  dynamics  of  rotating  elastic  structures  is  summarized,  showing  a  complex 
dependency  of  asymptotic  equilibrium  states  on  physical  parameters  such  as  elasticity  of  the  material  and 
total  momentum  of  the  structure.  For  multiply  articulated  rigid  bodies  undergoing  large  angle  rotations,  it 
may  also  be  shown  explicitly  how  there  can  be  a  sensitive  dependence  of  equilibrium  states  on  parameters 
of  inertia  and  angular  momentum.  This  suggests  that  recent  work  on  geometrically  nonlinear  beam  theories 
may  be  significant  not  only  for  transient  analysis,  but  for  accurately  modeling  asymptotic  dynamics  as  well. 


1.  INTRODUCTION 


Several  interesting  results  have  been  reported  in  the  recent  literature  dealing  with  the  nonlinear  mechanics  of 
rotating  elastic  structures.  Baillieul  and  Levi  (1987)  have  given  a  detailed  description  (by  means  of  a  geometric 
Lagrangian  formulation  of  the  basic  governing  equations)  of  the  asymptotic  equilibrium  behavior  of  several 
prototypical  structures.  For  one  such  structure,  consisting  of  a  rigid  body  with  a  simple  (Bernoulli-Euler)  elastic 

1  //v 

beam  attached,  it  has  been  shown  how  the  transient  dynamics  involve  an  interplay  between  inertial  forces  due 
to  rotation  and  viscoelastic  damping  of  vibrations  of  the  beam.  More  specifically,  it  has  been  shown  that  in  the  'or 
absence  of  external  forcing,  the  damping  will  have  a  dominant  influence  on  the  long  term  dynamics,  and  ultimately 
all  movement  of  the  appendage  with  respect  to  the  rigid  body  coordinate  frame  will  disappear.  Moreover,  in 


on- 


asymptotic  steady  state,  rotations  always  tend  to  a  constant  angular  velocity,  and  this  velocity  vector  is  aligned  - - 

with  a  principal  axis  of  the  steady  state  inertia  tensor  of  the  body- beam  complex.  ~j 

With  a  slightly  different  focus.  Simo  and  Vu-Quoc  (1987)  have  shown  that  in  modeling  the  dynamics  of  Codes 
_  and/or 
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rotating  beams  one  requires  the  use  of  geometrically  nonlinear  beam  theories  to  adequately  account  for  the 
influence  of  centrifugal  forces  on  bending  stiffness.  The  essential  feature  in  these  theories  is  their  incorporation 
of  higher  order  strain  energies.  While  the  geometric  formalism  and  the  simple  beam  model  used  by  Baillieul  and 
Levi  (1987)  incorporates  many  of  the  effects  of  inertial  forces,  it  does  not  capture  the  phenomenon  of  increased 
bending  stiffness,  and  this  suggests  that  such  a  model  must  be  refined  to  provide  a  faithful  representation  of 
transient  behavior  during  rapid,  large-angle  slewing  of  the  beam.  The  extent  to  which  more  complex  continuum 
mechanical  models  are  also  required  for  asymptotic  stability  analysis  is  the  subject  of  current  research. 

In  Section  2  of  this  report,  we  summarize  our  recent  work  on  the  asymptotic  equilibrium  dynamics  of  a 
rotating  structure  consisting  of  a  rigid  body  with  a  simple  elastic  beam  attached.  A  more  detailed  account  is 
given  in  Section  3  of  the  equilbrium  dynamics  of  a  kinematic  chain  being  rotated  about  an  axis  aligned  with  the 
force  of  gravity.  In  this  case  a  detailed  account  of  the  dependence  of  the  equilibria  on  the  various  mechanical 
parameters  points  to  a  high  degree  of  complexity  that  must  be  accounted  for  in  any  theory  of  the  mechanics  of 
rotating  structures  having  elastic  and  articulated  components. 


2.  EQUILIBRIA  IN  ROTATING  STRUCTURES 


In  Baillieul  and  Levi  (1987)  a  general  theory  of  rotational  mechanics  is  applied  to  a  detailed  study  of  the 
stability  and  equilibrium  dynamics  of  a  structure  consisting  of  a  rigid  body  with  a  cantilevered  beam  attachment. 
This  structure  is  depicted  in  Figure  2.1 


Figure  2.1:  A  rigid  body  unth  cantilevered  beam  attachment. 


Following  the  general  modeling  procedure  outlined  in  this  reference,  we  affix  a  ‘body  coordinate  frame'  so 
that  the  origin  coincides  with  the  point  of  attachment  of  the  beam,  and  the  zj-axis  is  aligned  with  the  axis 
of  the  undeflected  beam.  Elastic  deformations  of  the  beam  are  described  in  terms  of  a  function  ti(z.r)  which 
specifies  the  position  at  time  (  of  a  particle  whose  neutral  (or  undeflected)  position  is  (0, 0,  r)  in  the  chosen 
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body  coordinate  system. 

To  form  a  complete  model  of  the  dynamics  of  this  structure,  we  must  describe  the  time  evolution  of  the 
body  coordinate  system  with  respect  to  a  fixed  inertial  frame  together  with  motions  of  u(z,f)  with  respect  to 
the  body  system.  The  position  and  orientation  of  the  body  frame  with  respect  to  the  inertial  frame  is  described 
at  each  time  t  by  a  pair  Y(t)  €  50(3)  (giving  the  orientation)  and  y{t)  e  B3  (locating  the  origin  of  the  body 
frame  with  repsect  to  the  inertial  frame).  If  the  position  of  a  point  is  prescribed  by  a  vector  u  in  the  body  frame 
and  by  a  vector  V  in  the  inertial  frame,  then  u  and  V  are  related  by  U  ~  Y(t)u  +  y{t).  The  corresponding 
velocities  are  given  by 

L\  =  V(fiu  +  u, )  +  y 


where  fl(t)  is  the  skew  symmetric  matrix 


of  angular  velocities  about  the  corresponding  body  axes. 

With  these  notational  conventions,  we  may  prescribe  the  kinetic  energy  by 


T  =  iuT/u+mi,yI"yn?+  \mt\ly\ 

*  i 


l  f  4  u,)  t  y\  2dz. 

*  JO 


Here  u  =  (wj.wj.u'.)7  is  a  vector  whose  entries  are  the  above  angular  velocities,  c  is  the  center  of  mass  of  the 
rigid  body  in  the  body  frame,  mt.  is  the  mass  of  the  rigid  body,  we  have  scaled  the  linear  mass  density  of  the 
beam  to  be  one.  and  the  inertia  tensor  with  respect  to  the  body  frame  is  given  by 

/  “  lay  Jy  Iyg  . 

I*.  ly.  I. 

where  /,  is  the  moment  of  inertia  with  respect  to  the  x  axis  (=  zj-axis).  etc.  Since  the  beam  is  clamped  at 
the  origin  of  the  (2,,  z2.  zs)-coordinate  system  and  free  at  its  other  end.  the  following  boundary  conditions  are 
assumed 


u«(0. f)  =  -  (0.f)=  ~(/.f)  -0.  .  1.2. 
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and  i^a (0, / )  =  (l,t)  =  0.  These  boundary  conditions  are  standard  in  the  theory  of  clamped  free  beams 

Here  ut,  u2l  u3  are  the  deflections:  u  -  (uj,  u2, 2  -+  u3).  Note  that  u3  is  not  the  2-coordinate  of  u. 

The  equations  of  motion  for  our  rotating  structure  are  obtained  by  means  of  a  geometric  formalism,  the 
details  of  which  appear  in  Baillieul  and  Levi  (1987).  Thus  we  look  for  extremals  of  the  Lagrangian  L  =  T  -  V 
where  T  is  the  kinetic  energy  given  above  and  potential  V  is  an  appropriately  chosen  measure  of  strain  energy 
In  the  analysis  detailed  below,  we  focus  on  the  very  simple  strain  energy 

v(«)  =  l +  +  '  =  §-2, 

where  only  quadratic  terms  were  retained  and  the  material  is  assumed  to  obey  a  linear  Hooke's  law.  Here 
Hi  (respectively  /j2 )  gives  the  bending  elasticity  within  the  zz-plane  (yz-plane  respectively),  and  hi  is  the 
Hooke's  constant  giving  the  beam's  stretching  elasticity.  We  believe  the  asymptotic  analysis  based  on  this 
expression  gives  a  physically  realistic  picture,  but  Simo  and  Vu-Quoc  (1987)  have  pointed  out  that  certain 
important  features  of  the  nonlinear  dynamics  are  not  captured  by  so  simple  a  model. 

Viscous  damping  is  modeled  along  classical  lines  suggested  by  Kelvin  and  Voigt  by  introducing  the  dissipation 
function 

D  =  D(u.it)  =  ^  J\i(it'i)2  *  *2(«2)2  +  ks(it's)2di,  =  '  =  l,  (21) 

where  ti',',ti2  can  be  thought  of  as  the  rates  of  change  of  appropriate  curvatures,  while  ut  is  the  rate  of  change 
of  the  contraction  coefficient  u's.  The  k,  s  are  positive  constants  reflecting  the  rates  of  energy  dissipation  due 
to  deformation  of  material  in  the  beam. 

The  following  theorem  is  proved  in  Baillieul  and  Levi  (1987)  using  a  geometric  theory  of  Lagrangian  me 
chanics  with  damping. 

Theorem  2.1.  Given  the  tytlem  depicted  in  Figure  2.1  and  with  the  kinetic  energy,  potential  energy, 
and  dimpation  functions  described  above,  the  equations  of  motion  are  given  by 

Da  +  (nii,f  -t-  /  u)xy-,y=0  (2  2) 

Jo 

D3u  -t-  h^u  +  fcoi  •  Y~ly  =  0  (2  3) 
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where  the  quantities  in  these  equations  are  given  as  follows.  D(  )  =  !()  +  «*().  «(0  =  mo  * 
Jq  u  x  (£)u)<f2,  /  t«  t/i«  inertia  tensor  in  the  body  frame  defined  above,  m*,  t«  man  of  the  rigid  body 
component,  Y  —  VO.  and  d  is  the  differential  operator  defined  by  d  =  [d*,d*,  -3*),  p  —  diag(p \,P2,Ps) 
and  k  =  diag(ki,  k^.  ki). 

For  finite  dimensional  dissipative  mechanical  systems,  LaSalle's  invariance  principle  can  be  used  to  show 
that  states  asymptotically  approach  a  minimal  invariant  subset  of  the  zero  set  of  the  dissipation  function.  To 
some  extent  this  type  of  analysis  may  be  carried  out  for  dynamical  models  of  the  form  (2.2)-(2.4)  as  illustrated 
by  the  following  two  theorems 


Theorem  2.2:  Solutions  of  (2. £)•(£. J)  which  are  asymptotic  equilibria  (i.e.  solutions  which  also  satisfy 
D  =  0  with  D  as  defined  in  (2.1))  have  the  following  properties: 

(i)  there  is  no  dependence  on  the  time  variable  t  in  the  beam  function :  u(2,f)  =  tioc(')- 

(ii)  the  angular  velocity*:  is  a  constant  ; 

(Hi)  the  equilibrium  angular  momentum  is  a  constant.  ax  -  •/00u;00,  with  the  equilibrium  tnertia  tensor 
of  the  combined  body-beam  system  given  by 

./«  =  /+  ('  uTuE  -  uuTd;  -  m(ClCmE  -  CmCl)  (2  5) 

Jo 

where  E  -  the  identity  matrix  and  Cm  —  ^  (iTUc  -  /J  udz)  is  the  center  of  mass  of  the  body-beam  system 
(expressed  in  the  body  frame) 

(iv)  equilibrium  rotations  are  aligned  with  a  principal  axis  of  the  equilibrium  inertia  tensor,  and  thus 


Theorem  2.3:  The  asymptotic  equilibrium  beam  function  «/«(’)  and  the  asymptotic  equilibrium  angular 
velocity  wx  are  related  by  equations  (2  51  ),(2.6 )  together  with  the  fourth  order  system  of  ordinary  differential 
equations 

pdu=-nl(u  Cm),  (2.7) 

where  Cm  is  the  center  of  mou  of  the  body-beam  system  in  the  body  frame,  u  =  (t/j,«/2  “s  ■+  2),  and  Ylx  is 
the  skew  symmetric  matrix  whose  entries  correspond  in  the  appropriate  way  to  the  entries  in  wx  .  and  where 
the  boundary  conditions  are  as  prescribed  above. 
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(2.7)  prescribes  a  nonlinear  boundary  value  problem  since  depends  on  Uoo(-)  through  equations  (2.5) 
and  (2.6).  Some  idea  of  the  complexity  involved  in  explicitly  determining  ux(')  may  be  gleaned  from  our  analysis 
of  the  rotating  chain  system  in  the  following  section. 

3.  CLASSIFICATION  OF  EQUILIBRIA  IN  A  ROTATING  CHAIN 

While  a  primary  goal  of  much  of  the  work  reported  in  the  recent  liter¬ 
ature  on  the  rotational  dynamics  has  dealt  with  stability  issues  for  systems 
undergoing  three  degree-of-freedom  rotations  in  space,  valuable  insight  may 
also  be  gained  from  selected  restricted  problems  involving  single  degree-of 
freedom  rotations.  One  such  problem  involves  the  rotating  simple  kinematic 
chain  depicted  in  Figure  3  1.  Here  a  certain  number,  n.  of  links  of  various 
lengths  are  connected  to  form  a  single-strand  planar  kinematic  chain  We 
suppose  this  chain  is  suspended  in  such  a  way  that  the  force  of  gravity  tends 
to  extend  the  chain  to  its  maximum  total  length.  We  further  assume  there 

is  free  rotation  (with  no  actuation  or  friction)  about  all  joints  in  the  planar 
m, 

chain,  but  that  torque  may  be  applied  about  the  vertical  axis  (which  passes 
through  all  links  when  the  linkage  is  in  the  fully  extended  configuration)  An 
alternative  description  can  be  given  in  terms  of  the  Denavit  Hartenberg  (1955) 
parameters: 

(a,,a,,di)  =  (^.0,0) 

(a,,  a,  .d,)  =  (0, 1, ,  0) ,  t'  =  2 . n-fl. 

There  is  free  rotation  about  the  axes  , . . . , z„  depicted  in  F  igure  3.1.  and  an 

angular  rate  w  about  the  rP  axis  can  be  prescribed.  We  suppose  the  links  are 

massless,  but  each  joint  in  the  planar  chain  has  mass  m,.  )  (i  -  2 . n). 

Figure  3  1  "  degr*  an<j  there  is  a  mass  m„  located  at  the  tip  of  the  final  link.  The  angles  v, 

freedom  planer  kinemetrc 

chom  in  a  gravitational  depicted  in  Figure  3.1  measure  deviation  of  each  link  from  its  neutral  vertical 

field  rotated  ebout  its 

neutral  axis  position.  We  also  let  0O  denote  the  angle  of  rotation  of  the  mechanism  about 

the  vertical  axis  *0.  (Vo  >s  measured  with  respect  to  an  arbitrarily  chosen 
reference.) 


<j)  =  f 0 

c_> 


4  mn(2S2  4  ...  4  -  mr,S„g  -  0.  (3.l) 

where,  as  usual.  «,  =  sin  if>,  and  c,  =  cos  \p,. 

Before  showing  how  the  bifurcations  and  stability  characteristics  of  solutions  to  (3.1)  depend  on  various 
parameters  in  the  special  case  that  n  =  2.  we  present  some  general  results  to  the  effect  that  the  number  of 
distinct  solutions  (for  arbitrary  n)  is  not  less  than  2"  and  not  greater  than  4".  with  the  exact  number  of  solutions 
depending  on  the  values  of  the  parameters  and  w. 


Theorem  3.1:  For  email  values  of  )<*.■  | ,  there  are  precisely  2"  distinct  solutions  to  the  system  (S.l).  These 
solutions  are  given  explicitly  by  { ( v'  j ,  Vn )  '■  V-’i  —  0  or  »}.  The  index  (=  number  of  negative  eigenvalues 
of  the  Hessian  (|^))  of  any  such  solution  is  the  number  of  vector  entries,  ip},  which  are  equal  to  rr. 

This  theorem  is  proved  in  Baiflieul  (1987).  Physical  intuition  suggests  that  there  will  be  more  than  2" 
equilibrium  solutions  as  the  parameter  w  is  increased  in  (3.1).  To  count  the  total  possible  number  of  these 
solutions,  we  transform  (3.1)  into  a  systems  of  algebraic  (polynomial)  equations  by  means  of  the  trigonometric 
substitutions  x,  =  cos  v, .  y,  =  sin  .  The  resulting  equations  together  with  the  auxiliary  relations  r,2  ■+  yf  - 1  =  0 
(i  =  define  a  system  of  2n  quadratic  hypersurfaces  in  Cn.  We  have  the  following  Bezout  type  theorem 

regarding  these  algebraic  equilibrium  equations. 

Theorem  3.2:  For  generic  choice  of  parameters  mj , . . . ,  m„:  t\, . . . ,  tn\w  there  are  4n  (some  possibly 
complex)  solutions  to  the  algebraic  equilibrium  equations. 

Again,  this  theorem  is  proved  in  BaiHieu!  (1987).  Because  there  is  a  one-to-one  correspondence  between 
solutions  to  (3.1)  and  real  solutions  to  the  corresponding  algebraic  equilibrium  equations.  Theorems  3.1  and  3.2 
taken  together  imply  the  result  anounced  above  that 

Corollary:  Given  any  values  of  the  parameters  there  are  at  least  2n  but  no 

more  than  4"  solutions  to  the  system  (S.l). 

A  precise  description  of  the  way  in  which  the  numbers  and  stability  characteristics  of  solutions  to  (3.1) 
depend  on  the  parameters  in  the  problem  will  be  given  in  the  case  n  =  2.  It  is  convenient  to  reformulate  the 
equations  by  letting  l  —  t j,  l2  =  rt,  and  a  =  m2/(mj  +  m2).  If  we  moreover  assume  g  has  been  normalized 
to  be  1.  we  obtain  from  (3.1) 

(l  sin  •+  arfsin  V’2)  cos  ^ju>2  -  sin  =  0 

(f  sin  «/<  1  ■+■  r£sin  ^2)  cos  -  sin  1(2  =  0.  ( S.2 ) 

The  following  theorem  gives  a  substantial  answer  to  the  question  of  parametric  dependence. 

Theorem  3.3:  Let  a  and  r  be  fixed  in  the  system  of  equations  (S.2),  and  suppose  r  >  1 . 

When  ar  >  1,  there  is  a  monotonicaUy  increasing  dependence  of  the  number  of  solutions  to  (S.2)  on  the 
parameter  |w|,  and  as  |u>  increases  from  0  to  00  there  are  (counting  multiplicities)  f  .6,8,10,  or  12  solutions 
When  a r  <  1,  there  are  (counting  multiplicities)  4,6,8,10,12  or  16  solutions  to  (S.2).  In  this  ease, 
however,  there  will  not  in  general  be  a  monotone  dependence  of  the  number  of  solutions  on  the  value  0}  ja.1  . 


Figure  3.2  depicts  the  parametric  dependence  of  solutions  to  (3.2)  on  the  parameters  a  and  |uj  when 


r  =  1.5;  /  has  been  normalized  to  be  1. 


Figure  3.2  Bifurcation  locus  or  eauations  (3  2) 
in  the  case  r  =  15  The  top  curve 
is  cusped  at  approx  ec.33. 


The  tables  in  Figure  3.3  show  the  number  of  solutions  to  (3.2)  with  each  possible  index  of  the  Hessian  of 
L.  These  indices  of  course  give  the  dimension  of  the  unstable  manifold  associated  with  each  equilibrium. 


Figure  3.3(a):  When  or  >  1  the  number  of  critical  poinU  of  L  (=  solutions  to  fS.S))  increases  with  |w|  from 
4  to  IS.  The  table  gives  the  number  of  critical  point*  of  each  index. 


Figure  3.3(b):  When  or  <  1  the  number  of  critical  points  of  L  also  depends  on  |wl.  This  is  no  longer  a 
monotonie  dependence,  however,  as  indicated  by  the  eusped  bifurcation  locus  illustrated  in  Fig  S.  The  table 
gives  the  number  of  critical  points  of  each  index  eooresponding  to  each  possible  number  of  solutions  to  (S.S). 


4.  CONCLUSION 


A  detailed  account  has  been  given  of  the  way  in  which  the  equilibria  in  a  rotating  kinematic  chain  depend  on 
a  set  of  parameters  (which  are  easily  related  to  the  physical  parameters  of  momentum  and  inertia  in  each  link  in 
the  chain).  The  intricacy  of  these  results  points  to  a  corresponding  degree  of  complexity  that  may  be  expected 
in  the  equilibrium  dynamics  of  rotating  structures  where  elasticity  or  fluid  dynamics  must  also  be  accounted  for. 
For  rotating  elastic  structures,  we  may  ask  whether  the  higher  order  strain  energies  of  the  form  advocated  by 
Simo  and  Vu  Quoc  (1987)  will  lead  to  equilibrium  dynamics  which  differ  qualitatively  from  the  rotating  linear 
beam  models  studied  by  Baillieul  and  Levi  (1987). 

Near  term  research  goals  include  extension  of  this  type  of  equilibrium  analysis  to  articulated  structures  in  a 
zero  gravity  environment.  Also,  dynamic  simulations  and  empirical  studies  will  be  carried  out  to  obtain  a  more 
detailed  understanding  of  stable  equilibria. 
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